Self-consistent 3D radiative MHD simulations of coronal rain formation and evolution

Context. Coronal rain consists of cool and dense plasma condensations formed in coronal loops as a result of thermal instability. Aims. Previous numerical simulations of thermal instability and coronal rain formation have relied on the practice of artificially adding a coronal heating term to the energy equation. To reproduce large-scale characteristics of the corona, the use of more realistic coronal heating prescription is necessary. Methods. We analysed coronal rain formation and evolution in a three-dimensional radiative magnetohydrodynamic simulation spanning from convection zone to corona which is self-consistently heated by magnetic field braiding as a result of convective motions. Results. We investigate the spatial and temporal evolution of energy dissipation along coronal loops which become thermally unstable. Ohmic dissipation in the model leads to the heating events capable of inducing sufficient chromospheric evaporation into the loop to trigger thermal instability and condensation formation. The cooling of the thermally unstable plasma occurs on timescales that are comparable to the duration of the individual impulsive heating events. The impulsive heating has sufficient duration to trigger thermal instability in the loop but does not last long enough to lead to coronal rain limit cycles. We show that condensations can either survive and fall into the chromosphere or be destroyed by strong bursts of Joule heating associated with a magnetic reconnection events. In addition, we find that condensations can also form along open magnetic field lines. Conclusions. We modelled, for the first time, coronal rain formation in a self-consistent 3D radiative magnetohydrodynamic simulation, in which the heating occurs mainly through the braiding and subsequent Ohmic dissipation of the magnetic field. The heating is stratified enough and lasts for long enough along specific field lines to produce the necessary chromospheric evaporation that triggers thermal instability in the corona

The multi-thermal chromosphere: inversions of ALMA and IRIS data

Numerical simulations of the solar chromosphere predict a diverse thermal structure with both hot and cool regions. Observations of plage regions in particular feature broader and brighter chromospheric lines, which suggest that they are formed in hotter and denser conditions than in the quiet-Sun, but also implies a non-thermal component whose source is unclear. We revisit the problem of the stratification of temperature and microturbulence in plage now adding millimeter continuum observations provided by ALMA to inversions of near-ultraviolet IRIS spectra as a powerful new diagnostic to disentangle the two parameters. We fit cool chromospheric holes and track the fast evolution of compact mm brightenings in the plage region. We use the STiC non-LTE inversion code to simultaneously fit real ultraviolet and millimeter spectra in order to infer the thermodynamic parameters of the plasma. We confirm the anticipated constraining potential of ALMA in non-LTE inversions of the solar chromosphere. We find significant differences between the inversion results of IRIS data alone compared to the results of a combination with the mm data: the IRIS+ALMA inversions have increased contrast and temperature range, and tend to favor lower values of microturbulence in the chromosphere of plage. The average brightness temperature of the plage region at 1.25 mm is 8500 K, but the ALMA maps also show much cooler ($\sim3000$ K) and hotter ($\sim11\,000$ K) evolving features partially seen in other diagnostics. To explain the former, the inversions require the existence of localized, low temperature regions in the chromosphere where molecules such as CO could form. The hot features could sustain such high temperatures due to non-equilibrium hydrogen ionization effects in a shocked chromosphere - a scenario that is supported by low-frequency shock wave patterns found in the MgII lines probed by IRIS.Comment: 17 pages, 12 figures; accepted in A&A (added references, corrected typos

Accelerating Cosmic Microwave Background map-making procedure through preconditioning

Estimation of the sky signal from sequences of time ordered data is one of the key steps in Cosmic Microwave Background (CMB) data analysis, commonly referred to as the map-making problem. Some of the most popular and general methods proposed for this problem involve solving generalised least squares (GLS) equations with non-diagonal noise weights given by a block-diagonal matrix with Toeplitz blocks. In this work we study new map-making solvers potentially suitable for applications to the largest anticipated data sets. They are based on iterative conjugate gradient (CG) approaches enhanced with novel, parallel, two-level preconditioners. We apply the proposed solvers to examples of simulated non-polarised and polarised CMB observations, and a set of idealised scanning strategies with sky coverage ranging from nearly a full sky down to small sky patches. We discuss in detail their implementation for massively parallel computational platforms and their performance for a broad range of parameters characterising the simulated data sets. We find that our best new solver can outperform carefully-optimised standard solvers used today by a factor of as much as 5 in terms of the convergence rate and a factor of up to $4$ in terms of the time to solution, and to do so without significantly increasing the memory consumption and the volume of inter-processor communication. The performance of the new algorithms is also found to be more stable and robust, and less dependent on specific characteristics of the analysed data set. We therefore conclude that the proposed approaches are well suited to address successfully challenges posed by new and forthcoming CMB data sets.Comment: 19 pages // Final version submitted to A&

Prospects and challenges of numerical modeling of the Sun at millimeter wavelengths

The Atacama Large Millimeter/submillimeter Array (ALMA) offers new diagnostic possibilities that complement other commonly used diagnostics for the study of the Sun. In particular, ALMA’s ability to serve as an essentially linear thermometer of the chromospheric gas at unprecedented spatial resolution at millimeter wavelengths and future polarization measurements has great diagnostic potential. Solar ALMA observations are therefore expected to contribute significantly to answering long-standing questions about the structure, dynamics, and energy balance of the outer layers of the solar atmosphere. In this regard, current and future ALMA data are also important for constraining and further developing numerical models of the solar atmosphere, which in turn are often vital for the interpretation of observations. The latter is particularly important given the Sun’s highly intermittent and dynamic nature that involves a plethora of processes occurring over extended ranges in spatial and temporal scales. Realistic forward modeling of the Sun therefore requires time-dependent three-dimensional radiation magnetohydrodynamics that account for non-equilibrium effects and, typically as a separate step, detailed radiative transfer calculations, resulting in synthetic observables that can be compared to observations. Such artificial observations sometimes also account for instrumental and seeing effects, which, in addition to aiding the interpretation of observations, provide instructive tools for designing and optimizing ALMA’s solar observing modes. In the other direction, ALMA data in combination with other simultaneous observations enable the reconstruction of the solar atmospheric structure via data inversion techniques. This article highlights central aspects of the impact of ALMA for numerical modeling of the Sun and their potential and challenges, together with selected examples

Impacts of soil conditions and light availability on natural regeneration of Norway spruce Picea abies (L.) H. Karst. in low-elevation mountain forests

& Key message Natural regeneration of P. abies (L.) H. Karst. may reach high densities in lower mountain elevations. The highest densities were found in sites with moderate light availability, with low pH, and not near the riverbank. However, age-height classes differed in the predicted magnitude of response, but were consistent in response directions. Mosses and understory species typical of coniferous forests were positively correlated with regeneration density. & Context Norway spruce Picea abies (L.) H. Karst. in Central Europe is at risk under climate change scenarios, particularly in mountain regions. Little is known about the impact of environmental factors on the natural regeneration of P. abies in lowelevation mountain forests. & Aims We aimed to assess impacts of distance from the riverbank, soil pH, and light availability on natural P. abies regeneration. We hypothesized that (1) natural P. abiesregeneration would depend on light availability and soil pH and (2) there are understory plant species which may indicate the microsites suitable for natural regeneration of P. abies. & Methods The study was conducted in the Stołowe Mountains National Park (SW Poland, 600–800 m a.s.l.). We established 160 study plots (25 m2 ) for natural regeneration, light availability, soil pH, and understory vegetation assessment

Damping of coronal oscillations in self-consistent 3D radiative magnetohydrodynamics simulations of the solar atmosphere

Context. Oscillations are abundant in the solar corona. Coronal loop oscillations are typically studied using highly idealised models of magnetic flux tubes. In order to improve our understanding of coronal oscillations, it is necessary to consider the effect of a realistic magnetic field topology and the density structuring. Aims. We analyse the damping of coronal oscillations using a self-consistent 3D radiation-magnetohydrodynamics simulation of the solar atmosphere spanning from the convection zone into the corona, the associated oscillation dissipation and heating, and finally, the physical processes that cause the damping and dissipation. The simulated corona that forms in this model does not depend on any prior assumptions about the shape of the coronal loops. Methods. We analysed the evolution of a bundle of magnetic loops by tracing the magnetic field. Results. We find that the bundle of magnetic loops shows damped transverse oscillations in response to perturbations in two separate instances, with oscillation periods of 177 s and 191 s, velocity amplitudes of 10 km s−1 and 16 km s−1, and damping times of 176 s and 198 s. The coronal oscillations lead to the development of velocity shear in the simulated corona, which results in the formation of vortices seen in the velocity field that are caused by the Kelvin-Helmholtz instability. This contributes to the damping and dissipation of the transverse oscillations. Conclusions. The oscillation parameters and evolution we observed are in line with the values that are typically seen in observations of coronal loop oscillations. The dynamic evolution of the coronal loop bundle suggests that the models of monolithic and static coronal loops with constant lengths might need to be re-evaluated by relaxing the assumption of highly idealised wave guides

Origin of Norway spruce (Picea abies (L.) Karst.) in the Kartuzy Forest District

The origin of the Norway spruce in the Kartuzy Forest District was verified using the maternally inherited mitochondrial DNA marker mt15−D02 that shows distinct geographical distribution in Europe. We analyzed 462 trees from 26 populations. The majority of trees were of Alpine origin and they were present in all populations at a frequency of 6−90% (mean 59%). About 30% of trees were of the Carpathian origin and the least abundant (11%) were trees of north−eastern origin. The trees carrying north−eastern type of the mitochondial DNA marker were found only in 13 populations, but at a very different frequency (from 5 to 88%). The possibility of natural origin of Norway spruce in the area of Kartuzy Forest District was shortly discussed

Formation and heating of chromospheric fibrils in a radiation-MHD simulation

Aims. We examine the movements of mass elements within dense fibrils using passive tracer particles (corks) in order to understand the creation and destruction processes of fibrils. Methods. Simulated fibrils were selected at times when they were visible in a Hα image proxy. The corks were selected within fibril Hα formation regions. From this set, we selected a cork and constructed the field line passing through it. Other fibrilar corks close to this field line were also selected and pathlines were constructed, revealing the locations of the mass elements forwards and backwards in time. Finally, we analysed the forces acting on these mass elements. Results. The main process of fibrilar loading in the simulation is different to the mass loading scenario in which waves steepen into shocks and push material upwards along the field lines from locations near their footpoints. The twisted, low-lying field lines were destabilised and then they untwisted, lifting the material trapped above their apexes via the Lorentz force. Subsequently, the majority of the mass drained down the field lines towards one or both footpoints under the influence of gravity. Material with large horizontal velocities could also be elevated in rising field lines, creating somewhat parabolic motions, but the material was not generally moving upward along a stationary magnetic field line during loading. Conclusions. The processes observed in the simulation are additional scenarios that are plausible. The criteria for observing such events are described in this work. We note that it is desirable for our simulations to also be able to form more densely packed fibrils from material fed from the base of field footpoints. The experimental parameters required to achieve this are also discussed in this paper